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Noise poses a cliallenge for any real-world implementation in quantum information science. The 
theory of quantum error correction deals with this problem via methods to encode and recover 
quantum information in a way that is resilient against that noise. Unitarily correctable codes are 
an error correction technique wherein a single unitary recovery operation is applied without the 
need for an ancilla Hilbert space. Here, we present the first optical implementation of a non- 
trivial unitarily correctable code for a noisy quantum channel with no decoherence-free subspaces 
or noiseless subsystems. We show that recovery of our initial states is achieved with high fidelity 
(> 0.97), quantitatively proving the efficacy of this unitarily correctable code. 

PACS numbers: 42.50.Ex; 03.67.Pp; 03.67.Hk 



I. INTRODUCTION 

In the realm of quantum information science one of the 
most pervasive obstacles that must be overcome is the ef- 
fect of interactions with the environment - the resulting 
noise must be dealt with in any quantum system of prac- 
tical interest [l|, 0] . The theory of quantum error correc- 
tion has developed with this primary motivation in mind. 
The well-known passive error correction (or error avoid- 
ance) notions of decoherence-free subspaces (DFS) and 
noiseless subsystems (NS) 0, i, i, H, 0, i, Q have been 
extensively explored. Experiments have demonstrated 
the principle of DFSs in optical systems [13, HH, , ion 
traps [I^l , and NMR . The primary appeal of passive 
quantum error correction schemes is the fact they are 
passive - no correction operation is required once the 
noise acts. However, this occurs at the cost of a limited 
number of possible correctable codes for a given channel 
by relying heavily on symmetry in the noise jl5 |. 

When DFS/NS are not available or ideal, active quan- 
tum error correction (where a correction operation other 
than the identity map is needed) must be employed, re- 
sulting in an expanded set of correctable codes. Ap- 
proaches such as the stabilizer formalism [l^ provide for- 
mulas for finding and implementing codes in this regime. 
However, there can be a range of costs here that may neg- 
atively impact practical implementations. Clearly then, 
there are benefits to be had in the identification of cor- 
rection schemes that require the application of a recovery 
operation, but retain some of the appealing properties of 
passive codes. 

Here we experimentally implement a non-trivial exam- 
ple of such a correction scheme by applying random, an- 
ticorrelated noise to a two-qubit optical system. This 
noisy channel leads to the loss of encoded information 
and has no passive codes. However, we show that one 
qubit can be encoded against the noise such that the en- 
coded states are recovered with high fidelity by applying 
a single unitary recovery operation after each application 



of the noisy channel. This form of unitarily correctable 
code was recently introduced [13, [H, [3 ■ 



II. UNITARILY CORRECTABLE CODES 

One of the most general theoretical descriptions of how 
to deal with environmental interactions is given by the 
formalism of operator quantum error correction 17, 2(3]. 
In this theory, encoding, recovery and decoding of infor- 
mation takes place on subsystems of the full system. The 
key notion in operator quantum error correction is formu- 
lated as follows: a subsystem oiU = {H^ ® n^)(BK. 
is correctable for a quantum channel £ if there exists a 
recovery operation TZ such that for every pair of density 
operators p"^ G H"^ and G Ti.^ there is some other 
density operator e such that 

TZoSoVABip^ (Sa'') ^Vab{p^ (^T^)- (1) 

Here Ti,^ is an ancillary subsystem, K. is an ancillary sub- 
space orthogonal to Ti.^ CS) , Vab is the projection su- 
peroperator of H onto TL^^TL^ and by a quantum chan- 
nel, we refer to a completely-positive trace preserving 
map I']. Simply stated, a subsystem Ti,^ is correctable 
if, when we restrict our attention to it and anything it 
interacts with {Ti^ ® 'H^), there is a recovery operation 
that returns Ti^ to its initial state, while doing anything 

to n^. 

Decoherence-free subspaces and noiseless subsystems 
are two special cases of Eq. ([T]) , namely when the recovery 
operation is the identity map, and when dimT-^^ = 1 or 
dimTi^ > 1 respectively. Active schemes such as the 
stabilizer formalism are examples of Eq. ([1]) in generality. 

We say that a subsystem H,^ is unitarily correctable for 
a quantum channel £ if there exists a unitary recovery 
operation lA such that for every pair of density operators 
e and e there is some other density 
operator G TL^ such that 



Uo£o Vab{p^ ® a^) = Vab{p^ 



(2) 
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where W is a unitary map such that U : p ^ UpW, 
with U a unitary operator on ® . In the gen- 
eral case of Eq. ([1]), every recovery operation TZ can be 
implemented via a unitary map on a (in general) larger 
Hilbert space [Ij. Thus from another perspective, uni- 
tarily correctable codes are precisely the codes for which 
an extended Hilbert space is not required in the recovery 
process. 

It may not be immediately apparent why this is ad- 
vantageous as it is simply a special case of Eq. ^ . The 
first point to recognize is that Eq. Q is a generalization 
of the DPS /NS case since the identity map is a unitary 
map - thus the number of correction options for a given 
quantum channel will be at least as large. More impor- 
tantly perhaps, quantum channels with no passive codes 
can gain nearZj/-passive correction options. Secondly, this 
potential increase in correction options comes only at the 
cost of applying a single unitary map after the noise and 
no measurements. Thanks to these two points, unitar- 
ily correctable codes take advantage of two of the most 
appealing aspects of both active and passive correction 
schemes - increased range of correction options and easily 
implemented correction. 

From a practical standpoint, a correction scheme is 
only useful if one can find codes for which the scheme is 
capable of correcting a class of error operations of inter- 
est. In the case of unitarily correctable codes this can be 
achieved for the class of unital operations (those for which 
the identity operator is unaffected by the noise £"(/) — I) 
as shown by Theorem 2 of [l^ . The theorem states that 
the unitarily correctable codes for a quantum channel £ 
are precisely the DFS/NS for o 8, where ft denotes 
the dual map for f , defined via the relation between ex- 
pectation values: Tr(f(p)X) = Ty{p£^X)). O ne can 
readily check that £ is unital and trace preserving if and 
only if £^ is as well. 

The problem of finding the DPS /NS for a given chan- 
nel has been fully characterized [3, HH ■ In the case of 
a unital channel, the DFS/NS are the subsystems of the 
full Hilbert space that commute with each of the Kraus 
operators of the channel. This fact, together with The- 
orem 2 of [IS] provides a way to compute the unitarily 
correctable codes for a given unital channel. The next 
section includes an explicit calculation for a specific ex- 
ample. 



III. EXPERIMENTAL MODEL 

As a demonstrative example, we seek a simple quan- 
tum channel with no DFS/NS but with a unitarily cor- 
rectable code. To that end, consider the two-qubit phase- 
flip channel with Kraus operators ^^Zi, ^Z-^^ where 
Zi ^ az ® I and Z2 = I ® cFz- 

First note that there are no DFS/NS for this channel: 
f is a unital map so we must consider the operators that 



commute with both Zi and Z2. In the standard basis. 



Z^ = 



and, 



Zo = 
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Considering an arbitrary matrix in the Hilbert space. 



[ZuM] 



2mi,3 

2m2,3 

-2m3_i — 2m3,2 

-2m4,i — 2m4 2 



2toi,4 






(3) 



The off-diagonal blocks must thus be equal to zero for 
any element of the Hilbert space to commute with Zi 
and so all elements of the Hilbert space that commute 
with Zi must be of the form (for arbitrary a through h): 



f a b 

c d 

e / 

\0 g h 



Similarly, the elements that commute with Z2 are: 



f j k 

/ m 

n o 

\0 p q 



for arbitrary j through q. Therefore, the only operators 
that commute with both Zi and Z2 are of the form 



/ r 

s 

< 

VO M 



for arbitrary r through u. Since there are no non-zero 
off-diagonal elements we cannot encode any quantum in- 
formation - there are no DFS/NS for £. 

However, because £ is a unital map we can find the uni- 
tarily correctable codes for £ by looking at the DFS/NS 
for £^ o £. The Kraus operators for the unital channel 

£t o £ are ^■^I'^^j- In the standard basis 



Z1Z2 — 
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thus the operators that commute with the Kraus opera- 
tors oi £^ o £ are of the form 



/ a & 

c d 

e / 

\g h 



for arbitrary a through h. So, there are two one- 
qubit unitarily correctable codes for £ (the two one- 
qubit decoherence-free subspaces ior £^ o £): Ci — 
span(|00), 111)) and C2 = span(|01), |10)). 

Finally, in order to find a suitable correction operation, 
consider the effect of £ on an arbitrary two-qubit density 
matrix p = (pij): 



'-1 



|01)(01| + |10)(10| and also has dimension 2. Thus 



£ip) = 











-Pl,4 





P2,2 










-P3,2 


P3,3 













P4,4 



(4) 



There are many candidate unitary correction operations 
- for Ci the controUed-phase gate is one - it is also easy 
to see that both Zi and Z2 will correct either of Ci or 
C2 when taken as the unitary matrix of Eq. ([2])- In any 
experimental setting, some gates are easier to implement 
than others (e.g. in the optical setting, 2-qubit gates 
are difhcult because of the lack of photon-photon inter- 
actions). The flexibility in the choice of the unitary cor- 
rection operation is another appealing aspect of unitarily 
correctable codes. 

We have found what we set out to look for, a quantum 
channel with no DFS/NS, but two unitarily correctable 
codes. Although one qubit cannot be encoded with no 
recovery operation, one can be encoded with a single uni- 
tary recovery. Therefore this model allows for a clear ex- 
perimental demonstration of one of the key advantages 
of the unitarily correctable code approach ~ we gain a 
nearly-passive correction scheme for a quantum channel 
with no passive codes. 

An additional advantage of this model is that it can 
be realized in an optical setting. Defining \H) = |0) and 
\V) = |1), a phase-flip can be achieved by placing a half- 
wave plate (HWP) with its optic axis oriented along \H) 
in the path of a photon. Thus Zi {Z2) is implemented by 
a HWP in the first (second) photon path with nothing 
in the second (first). Since both Zi and Z2 will correct 
either codespace, we can experimentally implement the 
correction by placing a HWP in either photon path. 

Thus, the particular optical implementation of the uni- 
tarily correctable code we consider is as follows: consid- 
ering 



Uo£o Vab{p^ ® (T^) - Vab{p^ 



the full Hilbert space H. is the set of polarization states 
of two entangled photons. Ti. can be decomposed as H — 
{H^ ®n^)®JC where dimU^ = 1. The single encoded 
qubit will be supported on Pci = |00)(00| -f |11)(11|, so 
dimTi"^ = 2, and the ancilla space K, is supported on 



the projection superoperator Vab focuses our attention 
onto the codespace Ci : Vab '■ P ^ Pci pPci ■ The noise is 
provided by randomly fired anti-correlated phase-flips in 
both photon paths and has the form £ : p 1— > ^ZipZi -t- 
\Z2pZ2. Finally, the correction operation is provided by 
a single HWP placed in the second photon path, set to 
apply Z2, soU : p ^ Z2pZ2- 



IV. EXPERIMENTAL SETUP 
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FIG. 1: (Color online) Experimental setup, (a) Source of 
polarization-entangled photon pairs. An ultraviolet diode 
laser is used to pump a pair of orthogonally-oriented BiBO 
crystals to produce entangled photon pairs by spontaneous 
parametric down-conversion. The precise state produced is 
determined by the polarization of the laser light and can be set 
using a half- wave plate (HWP). A pair of compensation crys- 
tals are used to compensate group-velocity mismatch in the 
down-conversion crystals. The light is collected into single- 
mode fibers. More details can be found in the text, (b) 
Experimental implementation of the noise model, correction, 
and tomography. A quarter- wave plate is used to adjust the 
phase of the entangled pairs; together with the HWP in the 
UV laser and the fiber-based polarization controllers, these 
operations are sufficient to prepare an arbitrary pure state 
in C\. When the noise is on, the photon pairs are subject 
to computer controlled anti-correlated phase errors through 
liquid-crystal variable phase retarders (LCVR); a decision as 
to which LCVR will fire is made by a pseudo-random number 
generator every 1 s. The state of the light is measured us- 
ing quantum state tomography with a tomographically over- 
complete set of measurements and the maximum likelihood 
procedure [2^ . 
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FIG. 2: (Color online) Example set of experimentally mea- 
sured two-photon density matrices (left: real parts, right: 
imaginary parts), (a) The initial state was chosen to have 
a non-trivial value of both angles in Eq. ([Sjl ; the nearest pure 
state in Ci has angles 6 — 35.5° and — 46.5°. (b) The noise 
affected state has its fidelity to the initial state reduced to 
0.62; one can see that the primary effect of the noise is to 
flip the sign of the real and imaginary parts of the coherences 
rather than decohere the state, (c) The final state after both 
noise and correction. The fidelity with the initial state has 
been restored to 0.97 thus demonstrating the efficacy of the 
unitarily correctable code. 



Our experimental source is shown in Fig.fTJa). We use 
a 185 mW free-running UV diode laser (Newport model: 
LQC405-180E, centre wavelength 404.4 nm, bandwidth 
0.8 nm FWHM). To improve mode quality, the beam 
is passed through a spatial filter to obtain a 55 mW 
near-TEMoo beam. A motorized half-wave plate then 
rotates the pump polarization from horizontal to an ar- 
bitrary linear polarization. The diode laser pumps a pair 
of 0.5 mm thick orthogonally-oriented BiBO nonlinear 
crystals cut for type-I non-coUinear degenerate down- 
conversion with an 3° half-opening angle outside the crys- 
tal 23]. A 1 mm a-BBO and a 1 mm quartz both cut for 
maximum birefringence are used to compensate group- 
velocity mismatch in the down-conversion crystals. En- 
tangled photon pairs emitted by the source pass through 
bandpass filters (centre wavelength 810 nm, bandwidth 
5 nm FWHM) and are coupled into single-mode optical 
fibers. 



The light is coupled back into free space before appli- 
cation of the noise, as shown in Fig. [IJb). The noise £ 
is implemented by a pair of computer controlled liquid- 
crystal variable phase retarders (Meadowlark LVC-lOO, 
LCVR), one in each photon path. The LCVRs are both 
set to implement phase-flips. To simulate a noisy channel 
either one or the other LCVR is randomly fired with a 
switching rate of 1 Hz providing anti-correlated noise. 
The random switching was implemented using a soft- 
ware pseudo-random number generator (Lab View, Na- 
tional Instruments). Recovery is achieved with a half- 
wave plate oriented along \H) in the path of photon 2. 

We characterize all of our states using quantum state 
tomography. Fig [IJb) . There is a polarization analyzer 
in each arm comprised of a half-wave plate, a quarter- 
wave plate, and a polarizing beam-splitter. We use 
the tomographically-overcomplete polarization measure- 
ments {H, V, D, A, R, L} for each photon {R and L rep- 
resent right- and left-circular polarizations), i.e., 36 mea- 
surement settings on the pair of photons. Density matri- 
ces are reconstructed from the experimentally measured 
counts using the maximum- likelihood method (2^ . We 
performed quantum state tomography at each of three 
stages in the experiment: in the absence of noise or cor- 
rection, in the presence of noise, and in the presence of 
both noise and correction. 

Typical rates for the source are approximately 12000 
coincidence counts/s and 60000 singles counts/s when the 
fibers are directly connected to detectors. Using the half- 
wave plate placed before the down-conversion crystals, 
the fiber-based polarization controllers, and the tilted 
quarter- wave plate in the path of photon 1, the source 
can produce output states in Ci of the form 

^ cos2e\HH) +sm2ee"*'\VV). (5) 

With the half-wave plate oriented at 6* = 22.5° and the 
phase angle at (p — the system ideally produces the Bell 
state =^ + 1^^))- this case, we measure 

visibilities of 99.2% in the horizontal/vertical (H/V) ba- 
sis and 95.3% in the diagonal/ antidiagonal (D/A) basis. 
Reconstructing the state using quantum state tomogra- 
phy we find, in this case, that the fidelity 0, [2^ with 
0+) is 0.97, the tangle [20| is r = 0.9, and the hnear 
entropy [27| is Sl = 0.065. Thus our source is capable of 
producing highly entangled and nearly pure states. 



V. RESULTS &z DISCUSSION 

An example set of reconstructed density matrices, cho- 
sen to show the effect of the noise and correction on a 
state with unequal populations, and both real and imag- 
inary coherences, is shown in Fig. [51 Data was accumu- 
lated for 5 s per measurement setting. The initial state 
is shown in Fig. [DJa). This state has fidelity 0.98 with a 
state in Ci [Eq. (O] with 9 = 35.5° and cj) = 46.5°. Sub- 
jecting the photon pairs to the noise changes the state 
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product states \HH) and \VV), or HWP settings 0° (or 
90°) and 45°, respectively, are less affected by the noise; 
in the case of \HH) and |V^1^), the states are invariant 
under the noise. Those states closer to the maximally 
entangled states and 10") = -^{\HH) - \VV)), 

or HWP settings 22.5° and 67.5°, respectively, are most 
affected by the noise as it drives them to nearly orthog- 
onal states. The states are restored to greater than 0.98 
fidelity in all measurements thereby demonstrating the 
effectiveness of this unitarily correctable code across the 
entire codespace Ci. 



Qi w \y 

22.5 45 67.5 90 

Half-wave plate angle (deg) 

FIG. 3: (Color online) Measured fidelity between initial and 
final states. The mixed state fidelity [SJ was calculated on ex- 
perimentally reconstructed density matrices between the ini- 
tial states and noise-affected states (open circles, solid line 
theory) and the initial states and states affected by noise and 
unitary correction (closed, red circles). The data show that 
the effect of the noise depends strongly on the initial state. 
Initial states closest to the product states \HH) and \VV) 
are least affected while those close to the maximally entan- 
gled state \4>'^) are driven to nearly orthogonal states. The 
correction restores the photons to states with fidelity of higher 
than 0.98 with the initial states in all cases. 



to the one shown in Fig. [^I^b). The fidelity of this state 
to the initial state in Fig.[2lja) has been reduced to 0.62. 
Qualitatively, one can see that the effect of the noise 8 
does not decohere these states, but rather flips the sign 
of both the real and imaginary coherences. The correc- 
tion, hi, is implemented by placing a half-wave plate in 
the path of photon 2. When the noisy state was cor- 
rected, we obtained the state shown in Fig. ^c) where 
the fidelity with the initial state has been restored to 
0.97. The recovery of high fidelity states demonstrates 
the effectiveness of this unitarily correctable code. 

We characterized the effect of the noise and noise plus 
correction on a range of input states. We adjusted the 
source to states in Ci with real coefficients (0 = 0) and 
tuned those coefficients by adjusting the angle of the 
HWP in the pump laser, 0. The fidelity measurements 
of the noise affected states are shown in Fig. [3] where the 
initial states are shown as open circles and the theoret- 
ical prediction, F — 008^46*, is a solid black line. These 
data clearly show that the noise affects the states, but by 
different amounts. As expected, those states closer to the 



VI. CONCLUSION 

We have described a simple noise model which does 
not have a decoherence-free subspace or noiseless sub- 
system but does contain two unitarily correctable codes. 
We have implemented this noise model by constructing a 
quantum channel in a physical system, namely the polar- 
ization of a pair of optical photons, in which two liquid- 
crystal variable phase retarders fire in an anticorrelated, 
but random manner. We have sent photon pairs through 
this noisy channel in states ranging from separable to 
nearly maximally entangled and compared the impact of 
the noise and the noise plus unitary correction on the 
quality of the states. Our data shows that the noise 
can dramatically impact the fidelity of the output state 
with the initial, especially in the case of highly entangled 
states in Ci. However, the unitary correction restores the 
state quality with fidelities greater than 0.97 for all states 
in Ci. 

Our experiment shows how to translate the theory of 
unitarily correctable codes, a new method combining as- 
pects of passive and active quantum error correction, into 
experimental realization. As these codes expand the class 
of nearly-passive correctable codes, it is an interesting 
open question as to how these approaches could be uti- 
lized in more natural physical systems, as opposed to 
controlled application of noise. 
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